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We discuss the origin of decoherence in Josephson junction qubits. We find that two level systems 
in the surrounding insulator cannot be the dominant source of noise in small qubits. We argue that 
electron traps in the Josephson barrier with large Coulomb repulsion would give noise that agrees 
both in magnitude and in temperature dependence with experimental data. 
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Despite remarkable experimental breakthroughs of re- 
cent years QSIESIMIESIIIj the reahzation of quantum 
computer based on Josephson elements remains far away 
because developed circuits do not satisfy the stringent 
requirements on the amount of allowed decoherence im- 
posed by the quantum computation. Thus, it is crucial to 
identify (and possibly eliminate) the microscopic origin 
of the noise in these devices. 

It is commonly believed that the most important 
source of noise, at least for smaller qubits, is the fluc- 
tuations of the electric field generated by the environ- 
ment that are usually described in terms of the effective 
charge induced on a superconducting island. Recent ex- 
periments performed on single electron transistor (SET) 
and superconducting qubits provide an almost complete 
characterization of the charge noise power spectrum as 
follows. 

At low frequency (/ < lO'^ Hz) the charge noise spec- 
trum of Al — AlOx — Al SET shows a 1// dependence: 
Sq{uj) = ^ with intensity a = - IQ-^'e 
This level of the low frequency noise agrees well with 
the observed dephasiiig (free induction decay) in super- 
conducting qubits 0, S ■ Because its effect is greatly 
reduced by echo techniques |^, it is commonly believed 
that 1// noise does not extend above ~ 1 MHz. 
Recently, the direct measurements 0, of the tem- 
perature dependence of the l//-coefficient have shown 
the T^— dependence observed previously in different se- 
tups [11 [13. At high frequencies (/ > 1 GHz) the 
spectrum and the decay rate of the excited state of the 
superconducting qubits indicate the presence of a few 
quantum two level systems (TLSs) strongly coupled to 
the qubit 0, 0| superimposed on the Ohmic power spec- 
trum [1^ . Noteworthy, the crossover frequency between 
the high frequency linear noise and the interpolated low 
frequency 1// noise is of order of T, which suggests the 
same physical origin of the low and high frequency noise 

m 

In this Letter we discuss two microscopic mechanisms 
that might be responsible for the charge noise. First, 
we consider the effect of weakly interacting TLSs hosted 



in the amorphous material surrounding superconducting 
islands. We prove that this environment provides a sig- 
nificant source of decoherence but the detailed character- 
istics of its noise power spectrum are in qualitative and 
quantitative disagreement with the data; in particular, 
the observed dependence of o? requires unphysical 
density of states of TLSs Second, we discuss the 

effect of Kondo-like traps located inside the insulating 
Josephson contact. We show that these traps behave 
similarly to weak Andreev fluctuators. ^^l- We argue 
that for realistic parameters these traps produce noise 
very similar to the one observed experimentally. 

Quantum TLSs model. It is well established that TLSs 
dominate the low temperature behavior of most amor- 
phous materials at low temperatures [l^ Eol I2H . One 
usually assumes that these TLSs correspond to atoms 
that tunnel between two positions. Although their ex- 
act microscopic origin is not known, it is well established 
that they are well described by a pseudospin character- 
ized by dipole moment p — p = eda^ and controlled 
by the Hamiltonian Hq = Eaz + tax- Here E is the 
energy difference between the two minima, t is the tun- 
neling between them and those are distributed according 
to: 

P{E,t) = - — 10"3 (1) 

where v — IQ^'^ /cm^eV is the typical density of states of 
TLSs, d is the atomic distance and e denotes the dielectric 
constant of the medium. Different TLSs interact with 
each other via dipole-dipole interactions: 

H,nt = W (2) 

where is the distance between the i*'* and dipole. 
Notice that the effective strength of the interactions is 
given by the dimensionless parameter and is always 
very weak. This, for instance, implies that the effective 
field that other thermally active dipoles produce on a 
given one is small compared to the temperature. This 
small fluctuating field has, however, an important effect 
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on the coherent dynamics and is responsible for the low 
frequency 1// noise, as we show below. 

In order to estimate the coupling between a TLS and 
the qubit we notice that, in the charge basis of the qubit, 
its two states are characterized by the electrostatic poten- 
tial operator SVa^^i^ and thus produce electric fields that 
differ by 6£ = 5V/L at the position of the dipole, leading 
to the qubit-dipole interaction Hd-qub = S^P'^qub- ^^'^ 
characteristic length L is determined by the dipole posi- 
tion: it is of the order of the barrier thickness, 3 nm, 
for the dipoles in the barrier or L w 500 nm for the 
dipoles in the surrounding insulator. The resulting ef- 
fect on the qubit is conventionally described as the in- 
duced charge SQ = ne where the coupling constant 
K = p/eL ^ d/ L. Note that the coupling constant with 
the dipoles in the barrier is rather large, k ~ 10~^. 

Because the interaction term between the dipoles is 
small, it is convenient to work in the basis that diagonal- 
izes the Hamiltonian H^. In this basis Hq = ^ '^^j'^j 

The interaction between the 



where AE^ 



dipoles has two effects: it causes level broadening {I/T2) 
and it gives each level a finite decay time even in the ab- 
sence of phonons; we consider these effects in turn. We 
begin with level broadening. For this calculation we can 
neglect the back reaction and assume that dipoles fluc- 
tuate independently. In this approximation each dipole 
is subjected to the effective field h{t) generated by sur- 
rounding independent TLSs: 

Htll = E ''^ (^^^^ - (3) 

i 
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^ cosh^ (3AEj 



(4) 



where the coupling fc^ ~ — ' ■ Averaging 

ij 

Q over the distribution ^ we obtain the typical level 
broadening: 7 w (^) ^- expected, the dephasing 
level broadening is always much less than the tempera- 
ture. Note that this estimate assumes a significant to- 
tal number of thermally activated dipoles; this condition 
might not be satisfied in the thin Josephson barrier. 

The relaxation rate 7^ of dipoles is due to phonon 
emission or to the excitation of another dipole. The 
rate of the former decreases very rapidly at low ener- 
gies 7^ ^ AE^/ujjj where Ad ~ 10^ and becomes very 
low (~ 1 ms) for AE < 50 mK. The latter process dom- 
inates at very low temperatures. The time average of the 
relaxation rate of a given dipole can be found from Fermi 
Golden Rule: 
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HE' 



[AE,- AEjY +-f 



2 1 hj I ^ sin^ di sin^ 6j (5) 



This neglect of the higher order terms in the interaction 
can be justified if the dipolc-dipole interaction is weak 



compared to the line width {\kij\ < 7): more strongly 
coupled dipole pairs form new effective dipole modes and 
thus have to be excluded from the sum which is indicated 
by prime. Averaging jSJl we get for thermally excited 



dipoles: 7I « sin^ 0i i^^j T [23. 

As explained above, the effect of the fluctuating dipoles 
on the qubit is characterized by the induced charge noise 
spectrum: 



Sg{uj)= dt{5Q{t)5Q{Q))e"^' ^ {Key G{u:). (6) 

J —00 

with correlator G{uj) = ^. dt{al(t)al{0))e''^* . 

At high frequency the correlator is dominated by the 
resonant dipoles giving 
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Giuj) = E sin^ 0, 



{u; - 2AE,Y + 72 



(7) 



By substituting ^ into 10) and by averaging over the 
probability distributions we get 

eV 



Sq{uj) 



vp- 
e 



22^2 



(8) 



where V is the volume of the amorphous insulator. 

We conclude that the fluctuating dipoles give white 
noise spectrum at high frequencies. To estimate the noise 
level we first check that the amount of amorphous sub- 
stance is sufficient to provide a large number of thermally 
activated TLSs. In real qubits, TLSs are located in the 
barrier or in the surrounding oxide. For the barrier we es- 
timate V « l{f , so that the dipole density of states is 
about 0.1/ A'. Clearly this is not enough to give a contin- 
uous spectrum and we must conclude that the barrier can 
contribute to the charge noise at high frequency by pro- 
viding only few TLSs, strongly coupled to the qubit; this 
might at most account for resonance peaks. The situation 
is different for the surrounding oxide which has volume 
V > 10^^"^, in this case the number of thermally acti- 
vated TLSs might be sufficient to produce a continuous 
spectrum of relaxation times with Sq{ui) k, 10^^e2 K^^ . 
Experimentally, the linear high frequency noise and the 
low frequency 1// noise both extrapolate to the same 



value 



10- 



at w - T « 120 mK M 



This 



value is very close to our estimate showing that the TLSs 
might provide a non-negligible source of noise at these in- 
termediate frequencies. 

At low frequencies the noise is dominated by thermally 
excited dipoles that behave as classical fiuctuators at fre- 
quencies less than their relaxation rates: 



E cos^ I 



(9) 



The relaxation rates 7} have a wide distribution, 
P(7) 1/7, producing upon averaging, 1// noise: 
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Its linear temperature dependence is a direct consequence 
of the constant density of states that results in a linear 
dependence of the number of thermally activate TLSs: 
n-TLS = wPo^ where W w leV is the energy scale of 
the disorder and po = lO^^/cm"^ is the bare density. Sim- 
ilarly to the high frequency noise, this noise is of the 
right order of magnitude at lowest temperatures but 
grows linearly with temperature instead of following the 
dependence. Moreover, the phenomenological as- 
sumption that the number of thermally active dipoles 
is Utls = ^nPo^ h^^s an other problem: the observed 
value of 1// noise at low frequency implies that the total 
density po = ^Po ~ 10® po is much larger than one per 
atom. Thus, in order to account for the data, one needs 
a physical mechanism characterized by a small energy 
scale W. 

Kondo-like traps model. A natural candidate for the 
small energy scale is the superconductive gap A. In the 
bulk of superconductor all electron states are gapped and 
low energy electron traps seem unlikely. It is very likely, 
however, that numerous electron traps appear in the in- 
sulator at the SI boundary, with energies distributed in a 
broad energy range. The Coulomb repulsion between two 
electrons in the same trap is very strong prohibiting dou- 
ble occupancy. In the absence of such repulsion, the ef- 
fect of the electrons hopping from one trap to another and 
their hopping via Andreev process would give the correct 
frequency and temperature dependence [l^l , but requires 
an unphysical density of states of these traps, similar to 
the TLS discussed above. The strong on-site repulsion 
leads to a formation of a weak Kondo-like resonance at 
the Fermi surface and provides the desired low energy 
scale. This physics is described hy H = Hd + Hgd, where 
Hd and Hgd describe the Hamiltonian of the Anderson- 
like traps and of the coupling between the quantum im- 
purities and the superconductor: 

■ia i 
Hsd = "^"^ hkicl^Cdicr + h.C.] , (11) 
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'^dicr 



^dit 



i'^dia) creation oper- 



Here ndia 

ator of a conduction (localized) electron and U is the 
on-site repulsion energy. In the absence of supercon- 
ductivity, the Kondo effect jl^ leads to the appear- 
ance of the resonance at the Fermi energy with width 



where F, = t" 



VAl 



and VAl is the density of states in the Al. This reso- 
nance is only weakly affected by the superconductivity 
if Tk ^ A. In this case, the ground state of the system 
(impurity plus superconductor) is a singlet and contains 
even number of electrons, while the first excited state is 
localized with the energy e^f < A. In the opposite limit 
Tk <C a the trap is always occupied, the ground state is 
a doublet and contains odd number of electrons. Because 



the full Hamiltonian does not mix states with different 
parity, as a function of Tr-/A a singlet and doublet level 
should cross at k, A. Numerical studies confirm this 
and give the value « 0.3A. Close to this value we ex- 
pect that the energy of the excited state varies linearly: 



CK = ^ {Tk -T*k) T, 



K 



0.3A 



(12) 



with I? ~ 1. Thus, the electron traps characterized by 
Tk ^ T^ provide a density of localized states at very 
low energy. Since the Kondo temperature depends expo- 
nentially on the energy e° of the magnetic impurity, the 
density of these low energy states is 



(13) 



v{e) = poV / dp{eK)S{eK ~ e) 

JO 



Tk 



K 



where po is the bare density and V is the volume of the 
insulator barrier that contains traps. The density of these 
states is much larger than the one of the traps without 
on-site repulsions. However, their effect on the qubit is 
not so large because it is compensated by a small weight 
of the Kondo resonance w = Tk/s^ 

The charge of the qubit produces different electrostatic 
potentials on the traps, so electrons/holes tunneling be- 
tween them causes the qubit decoherence. We distinguish 
fast processes in which electrons move a short distance 
r < ^ which are responsible for high frequency {to > T) 
behavior of S{to), and exponentially slow ones where the 
distance is large r ^ £^ (^ is the coherence length of the 
superconductor). If the size of the contact, /, is small 
{I «C C), the relaxation might be dominated by tunneling 
from the traps inside the Josephson barrier to the one in 
the surrounding insulator that has much larger volume; 
however for a typical contact / ~ ^ we shall ignore the 
distinction between the traps in the barrier and in the 
insulator nearby. In order to estimate the contribution 
of the fast processes, we compare it with the one of TLSs 
There are two important differences: first, at a given 
frequency uj any pair of traps that differ by uj in energy 
can play the role of the two level system if one of these 
traps is empty and another is occupied; so the density 
of states of these pairs is u!i''^{e). Second, the effective 
dipole moment of the pair contains the additional small 
factor, w: p = wea where a is the typical distance sep- 
arating the trap from the superconductor. Combining 
these factors we get 



Sg{uj) = {W 



L 



2 ( P0V\ 

\W*k) 



e^uj 



(14) 



The electron tunneling amplitude falls off exponen- 
tially for traps separated by distances larger than ^ lead- 
ing to the logarithmic distribution of long relaxation 
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times: dP{T) = dr/r. Neglecting the logarithmic depen- 
dence of the insulator volume contributing to the slow 
dynamics we estimate the number of thermally active 
pairs by (Tz/(e))^. The dipole moment corresponding to 
this transition is wea, so we get S'q(w) ~ ^ with 



a w Kwe 



(15) 



As expected, the linear frequency dependence at high 
Lo translates into a T^-dependence of the low fre- 
quency power spectrum. Assuming, as before, that 

V w 10"^!^, po ~ lO-^i^ and w w 10"'' we estimate: 
a « 10-2 - IQ-'^e for T ~ 100 mK H^. We conclude 
that these fluctuators, although very weak, compensate 
their weakness by their large number, providing a value 
for the intensity of the 1// noise in reasonable agree- 
ment with the experiments. Furthermore, they give the 
correct dependence of the noise power at low and high 
frequency. Notice that the electron moving in and out of 
the trap has also a large effect on the Josephson current 
if this trap is located in the Josephson contact, so these 
fluctuators might be responsible also for the noise in the 
critical current fluctuations as observed in 01 ■ 

In order to estimate the upper cutoff. A, for the 1// 
noise we compute the tunneling amplitude between two 
traps separated by distance ^ located in a box of size ^. In 
the absence of Coulomb repulsion on each trap, the tran- 
sition amplitude between two traps with energy e reads: 
A = J2a g'l'^"*'"''**^ ■ Assuming that the wavefunctions 
of the electrons of the superconductor are random we find 
the tunneling amplitude |Ap = /i'aiS,^^- The on-site 
repulsion changes the energy, e, of each trap to H12|l and 
replaces F — giving 



\Ao\' = 



(16) 



where 6 = 1/(1^^;^"^) is the average distance between the 
levels in the box. For a typical Al island Aq « IG^^T^, 
which is smaller than experimental temperatures. For 
a typical aluminum island (A^i« 200fieV) we estimate 
Aq « 10® Hz. In fact these processes involve energy 
transfer to a thermal mode, so A < ^o- 

Conclusion. We have discussed two possible micro- 
scopic mechanisms responsible for the charge noise in 
qubits. We found that for realistic values of the param- 
eters, the conventional TLSs in the insulating Josephson 
barrier cannot provide a significant source of decoherence 
for small qubits, although they are likelyto be the dom- 
inant source of noise for large junctions [3|. The TLSs in 
surrounding insulator might provide a significant source 
of noise at uj <r ^ 100 mK, but would give frequency in- 
dependent Sq{uj) at large uj and linear T— dependence of 
the 1// noise. We considered the effect of Anderson-like 
traps in the Josephson contact and found that quasipar- 
ticle hopping between these traps provides noise of the 



right order of magnitude with the correct temperature 
and frequency dependence. In order to test this mecha- 
nism one needs to develop a more detailed analytical the- 
ory for the maximal rate of the thermally activated hop- 
ping between the traps and for the effect of the magnetic 
field that suppresses the superconductivity and compare 
these predictions with the data. 
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